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Discrete Hybrid Polyharmonic Cubature 
Formulas with weight on the disc, with error 
bounds 


O. Kounchev and H. Render 


Abstract 

This paper is a second part of our study of the Discrete Polyharmonic Cubature 
Formulas on the disc. It completes our study and provides a satisfactory cubature 
formula in terms of precision and number of evaluation points (coefficient of efficiency 
of the formula). 

1 Introduction 

The present paper is a continuation of paper |43j . where we have provided a Discrete 
Polyharmonic Cubature Formula for integration of functions multiplied by a weight on 
the disc. The main merit of this formula is that unlike majority of research in the area, we 
provide constructive error bounds. However an issue arises with the number of evaluation 
points. In the present paper we show how to solve this issue by using spline approximation 
methods. This approximation increases the error of the evaluation of the integral but it 
is again under control and we show that it is a reasonable compromise which gives much 
better result compared to the alternative methods. 

We decided to provide a separate presentation of the spline version of the Discrete 
Polyharmonic Cubature Formula, which we call here Hybrid Discrete Polyharmonic 
Cubature Formula, in order not to overburden the attention of the reader of the original 
paper |33|. We will rely upon the notations and definitions provided in [33]. 

As in [43] . we consider numerical evaluation of integrals of the type 

Iw(f)= [ f(x)w(x)dx (1) 

Jd r 

where w ( x ) is a (not necessarily non-negative) weight function on the disc Dr in the plane 
M 2 . As in |33] we provide experiments with the weight functions 

w (1) (x,y) = and w (2) (x,y) = \y\ . (2) 

V x + V 

We provide the experiments in a form facilitating comparison with the results in [33]. 

For completeness of notations in the present paper, we repeat the main formulas used 
in the previous paper [35] : 
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1. The Fourier series of w is given by 


oo a k 

W (re l(p ) =w(r cos ip, r sin ip) := Y^2 w (k/) ( r ) Y {k,l) O) ■ ( 3 ) 

k=oe=i 

2. Orthonormalized spherical harmonics are defined by 

*(o,i) M = l/V^F (4) 

Y (k,i) (<p) = cos kip and Y^) ( <p ) = T s i n fc<p. (5) 

for integers k > 1. Then Yj^^) is an orthonormal system for k > 0, £ = 1,<Zfc, 
where a*, = 2 for k > 1, and ao = 1. 

3. The (A;,£)-th Fourier coefficient of a complex-valued continuous function / (re 1 '?) is 
given by 

/*27T 

f(k,i) O’) ■= f (re lif ) Y (M) (ip) dip for k > 0, l = 1, a k (6) 

Jo 

4. The corresponding Fourier series of / is 

oo a k 

f (re 1 *) = f(r cos ip, r sin ip) := O') Y (M) (ip). (7) 

k=0£=1 

5. The Fourier series of a polynomial P (x) is of a very special form: there exist poly¬ 
nomials P(k,e) an( l a number N < deg P (x) such that 

N a k 

P(x) = P (r cos ip, r sin ip) = (n 2 ) r k Y (k/) (ip); (8) 

k=0£=1 

the representation ((HJ) is called Gauss decomposition or Almansi expansion of a 
polynomial p. 

6. Hence, the Fourier coefficient P( k ,£) ( r ) °f a polynomial P (x) is of the form 

P(k,£) ( r ) = P(k,£) {r 2 ) r k . (9) 

7. The integral ©> in polar coordinates, becomes 

p2tt nR 

I w (/) = / / f (r cos ip, r sin ip) ■ w (r cos ip, r sin <p) ■ r dr dip. 

Jo Jo 

8. We obtain 

p2tT pR 00 a k pR 

4- (/) = / / f (re lip ) w (re lip ) rdrdp = VV / f k ,£ (r) uj( k/) (r) rdr. (10) 

Jo Jo k=Q£=1 Jo 
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9. For the constant weight function w (x,y) = 1 = y/27rYr 0 u one obtains simply 

p2n pR _ pR 

I\ (/) = / / f (r cos tp,r simp) rdr dtp = V2 tt / /(o,i) ( r ) r dr. (11) 

io io JO 

10. We shall assume that each function iV(k,t) ( r ) does not change the sign over the 
interval [0, R] , a property which we call pseudo-definiteness 


W(k,i) ( r ) > 0 on [0, /?] or 

W{k,t) ( r ) <0 on [0, R] . 


( 12 ) 


11. By the Almansi formula ([5)1 . and a change of the variable p = r 2 , the one¬ 
dimensional integrals in (TTTT1) used for computing I w (p) are equal to 


pR pR 

/ Pkjt ( r ) W(k,l) ( r ) rdr = / Pk,£ {r 2 ) r k w {k)tj (r) rdr 
Jo Jo 


1 


rR 2 


P{p)-p k/2 W {k/ ) (y/p)dp. 

(13) 


12. Assuming [12] to the integral m with measure p kd2 W( k jo {y/p) , we apply the Ap¬ 
point Gauss-Jacobi quadrature: we obtain the nodes < ••• < *jv,(fe,£) an d the 

weights o,..., \N,(k,i) which are either all positive or all negative. Due to the 
exactness of the Gauss-Jacobi quadrature for any integer 0 < s < 2N — 1 we obtain 
the equalities 


N 


rR 2 


E A h(M) • fS Uk,e) = 2 P S P k/2w (k/) WP) dp 

j=1 Jo 


rR 


r 2s r k w (k/) (r) rdr 


(14) 


13. Hence, for a polynomial /, for which deg ,/A.n < 2N — 1 for all (k,£), by using the 
Gauss-Jacobi quadrature ([TO . the integral (fTUD becomes 


oo a k R 

Iw (/) = / fk,e (r) (r) rdr 

k=oe=i Jo 

OO a k R 2 


(15) 


1 ^ /*ri r . 

= 2^E / /*.< (V^) P"* /2 (p^Mm) (v/p)} ^ 

Z fc=0^=ld° 

= ^ ly (/) ■ 

where we have put 

oo a fc tv 

J N° ly (/) : = 5 EEE*i,(M> • Em) ' f(k,e) iV^Jk /)); (16) 


k=0Z= 1 J=1 


14. As the values fk/ are Fourier coefficients, they can be approximated by 

means of Discrete Fourier transform of the function /. Thus for every two integers 
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K, M > 0 we introduced the Discrete Polyharmonic Cubature with parame¬ 
ters (. N,M,K ) in the paper [[53], by putting 


K a k N 


( nJm,k) (/) 2 £££a, ( m) ’ Akm) ’ Am) (vAAm)) 


(17) 


7T 

M 


fc=0£=lj=l 
K a k N M 


’ A‘ : (M) ■ Am) f M 1 

fc=0£=lj=ls=l V ' 


/ ( \A?',(M) e 


27TS 

M 


(18) 


where the Discrete Fourier transform is given by 


f (M) 


(r) 



(19) 


15. Its coefficients (weights) | • t,- ■ Y(k,t) (i&r) f have varying signs but they 

satisfy the following remarkable inequality 


sEEEE 

k=oe=ij=is=i 


K a k N M 


f 2irs 

Aj,(M) ■ tj,(k,e) ’ Am) y-j^f 


< y/TT ||u?|| , 


( 20 ) 


where we have assumed that the weight w satisfies 

°o a k ft 


u> := 


^ r-ti 

-■= X] Yh / Km) ( r )\ rdr < °°- 

fc=0£=1 do 


( 21 ) 


This inequality is proved in |43| , by an application of the famous Chebyshev extremal 
property for the Gauss-Jacobi quadrature: 

JA _fc i rR 2 rR 

2A I A,(M) | A(M) < g J Q I W (M) (v^) I d P = y o |™(M) (A | rdr. (22) 


16. Let us note that formula (1181) needs (2 K — 1 )-N-M evaluation points. The coefficient 
of efficiency of a cubature formula may be introduced in the following way (see e.g. 
[66]): If n is the number of nodes of the quadrature, m is the number of linearly 
independent functions which are integrated exactly, and d is the dimension of the 
Euclidean space, then the coefficient E is defined as 


n (d + 1) 

For Gaussian quadratures in one dimension, E is obviously equal to 1. For the 
Discrete Polyharmonic Cubature formula (1171) we need precisely (2 K — 1) • N ■ M 
sampling points - hence 2 coordinates for every point and 1 coefficient. On the other 
hand, as we have seen in [43j, the subspace of exactness of formula (HZD has dimension 
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2N x (2 M — 3 — 2 K), hence the efficiency coefficient of the Discrete Polyharmonic 
Cubature formula is 

2TV x (2 M - 3 - 2 K) 

3 (2/i -1)xNxM 

which is much less than 1, hence is prettily bad. The main purpose of the present 
paper is to show how by using spline approximation and sacrificing some precision, it 
is possible to reduce the number of point evaluations which will improve essentially 
the coefficient E. by making it very close to 1. 

Now we come to the main point of the present paper. 

Since the point evaluations of a function / can be very costly, it is of advantage to 
reduce their number by sacrificing some precision. There are simple ways to decrease 
essentially the number of evaluation points, e.g. by applying spline approximation to the 
coefficients f(kj) (\Aj,(M)) • F° r example, let us assume that for some integer N\ > 1, 
we are given the radii Rj, j = l,2,...,Ni with 0 < R\ < R 2 < ■ ■ ■ < Rn^ < R, and the 
N\ x M functional values be given as well: 

/ (Rje ilfs ) for s = 1,2,...,M, and j = 1, 2,..., N x . 

Thus we will have only Ni x M sampling points of the function / which is the usual frame¬ 
work for the cubature formulas. Then, by formula (1191) . we will find the Discrete Fourier 
transform approximation fj.f'P ( Rj ) to j\ t f {Rj) , and we may use interpolation splines 

with data {fk,i (Rj)}f= 1 to approximate all values fjpp f° r 3 = b 2,..., N. The 

resulting formula introduced below is termed (Spline) Hybrid Discrete Polyharmonic 
Cubature formula (1241) . 

As with the experiments in [ 13] , the experiments with formula (1241) . show excellent 
performance. This is due to the fact that we are able to control the error bounds of 
formula (|24l) . where we now have the parameter N\, for the knot evaluations, see Theorem 
[3] below. 

The plan of the paper is the following: In Section [2] we define the spline based Hybrid 
Discrete Polyharmonic Cubature and provide error bounds for it. 

In Section 0 we provide experimental results for the discrete polyharmonic cubature 
with respect to the first weight function a/ 1 -* in (f2j) . 

Section H] contains experimental results for the second weight function w^ 2 \ 

2 The Hybrid Polyharmonic Cubature formula 

Let let us assume that on the two-dimensional disc Br of radius R, the values of the 
function / (x ), x € Br C K 2 are given on a regular grid: We assume that the integers 
N\ and M are fixed in advance, where M is odd. Then we consider the concentric circles 
with radii Rj, for j = 1 , 2 ,..., N\ , with 0 = Ro < Ri < Ro < ■ ■ ■ < Rn± = R, and we 
assume that the values of the function / are given: 

f (R m e l ^i s ' s j for m = 1,2,..., N x , s = l,2, 

Formula (fl9l) provides an approximation fpp (. R m ) to the functions ( R m ) • 
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The value fk/ (r) at the knot r = 0 is determined from the following arguments: In 
the case of k > 1, by Proposition 1 in [9], we know that if the function f € C°° (Dr) then 
f(k,i) € C°° [0, R] , and 

f(k,t) (0) = 0. 

Hence, we will put also 

/(M) (°) = 0 for k ^ L 

On the other hand, we see by formulas (HD, Q that the very first coefficient satisfies 

/(o, 1 ) (°) = ^/(°), 

hence, we will have 

In general, if V € C Nl+1 is a vector, we will denote by 


SPL [V] (r) (23) 

the value at r, 0 < r < 1, of the interpolation cubic spline with knots {Rj}Pl 0 an d i n t er P°- 
lation data which satisfies the not-a-knot boundary conditions, see Proposition 

[1] below. If we denote by S m (r) the not-a-knot interpolation spline which satisfies 


Sm (-^s) — &m,s 

(where 5 denotes the Kronecker symbol) then we obtain the representation for the data 
V, given by 

Ni 

SPL [V] (r) = S m (r) V m . 

m =0 


The idea is now to approximate the values fjpp through the not-a-knot 

spline SPL interpolating the values | fgPi) (hy) j at the nodes 0 = Rq < Ri < R? < ■ ■ 


3 = 0 

• < = R. Thus, from the Discrete Polyharmonic Cubature with parameters (IV, M, K ) 

defined in (ED , after substituting ( y/tjjk/)) with SPL j (R m ) |( y/tj,(k,e) 

we obtain the Hybrid Polyharmonic Cubature formula: 

K a k N 


7-spline 

1 (N,M,K,N 1 ) 


(/) 9 

w k=oe=ij=i 


x t. /I n, x SPL 
mM) 




Obviously, we have the equalities: 

K a k N Ni _ k _ 

I(N,M,K,Ni) (?) = 2 X/ X/ X ‘S'm (v^(M)) ^(M) ( 25 ) 

fc=0 f=l j=l m=0 

JVi M 

m=0 s=l 
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where we have put 


1 \ n ( / -\ 27T ( 2tt \ — ^ 

Wm,s '■= 2 Z_^ W t L(k/)) Tf Y {W) [ Jff S ) X 

k= 0 1= 1 i=l v 7 

The last formula (|25l) shows that I^'mk Ni) (/) * s a cubature formula in the usual sense, 
having IVi • M nodes. 

Let us note that formula (1251) has the same knots as the Midpoint cubature rule, see 


Let us put 

h := max (Ri+ 1 — Rf) • (26) 

i 

Then the error bound for not-a-knot spline is given in the following Proposition. 

Proposition 1 Let g € C 4 {[a,b\). Assume that the knots of the interpolation are a = 
Rq < R\ < ■ ■ ■ < Rn 1 = b. Let N\ > 5. Let s ( t ) be the not-a-knot intepolation cubic 
spline, i.e. the C 2 piece-wise cubic polynomial function which satisfies the interpolation 
conditions 

s{Rj) = g(Rj ) for j = 0 , 1 ,..., N\ 

and the not-a-knot conditions at the two penultimate knots: 

s (3) (Ri - 0) = s (3) (R 1 + 0) and s (3) (R Nl _ i - 0) = s (3) (R N i_i + 0). 


Then the error in the interpolation satisfies 


ig-s) 


(r) 


< C r 


(4) 


x h 


4—r 


for r = 0,1,2,3, 


^co fo] 

where the constant C r is independent of the function g and the size of the grid h. 


For the proof we refer to [12] and [10J, inequality (2.1). 


Remark 2 The proper value of the constant Co above may be estimated by the some 
arguments, provided in m, p. f22, and Theorem 4.4.8 therein. Thus the value as provided 
in Theorem 4.4.8 m seems to be a reasonable approximation: 

C 0 < —• 

Now we may prove the following main theorem. 


Theorem 3 Let the function f € C°° (-Dr) . Let N\ > 5. Then the difference 


jspline / « _ jpoly / n x 

1 (N,M,K,N 1 ) w > 1 (N,M,K) W > 

between the Hybrid Polyharmonic Cubature formula (24\ ) and the Discrete Polyharmonic 
Cubature formula m satisfies the following estimate: 


jspline 

1 (N,M,K,N 1 ) 


(/) - If* 


(. N,M,K) 


(/) < V^Ch 4 X 


< VnCh 4 X 


d 4 f ( re llf 
dr 4 

d 4 f (re* 
d r 4 


1 

) 


K a k R 

__ *££/ |™(M) W| r dr 

L x (Df) k=0£=i Jo 
x ||?n|| . 

Loo (D r ) 
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Here C > 0 is the constant Co, independent of the function f and the parameter h, as 
provided by Proposition QJ 


Proof. In order to find the error bound of this approximation we recall that for every 
r with 0 < r < R we have by definition m, 


/£?(') = ££/ 


2vr 


M 


M 


s=l 


re*® 8 ) Y, 


(M) 


2vr 

m‘ 


By the linearity of the spline interpolation, for every r with 0 < r < R we obtain 

M 

C an -i JVi Vtt 

SPL 


{C> [{/ K e ‘ Ss )} 

J S=1 L 


N! 
m =0 


2tt 

( r ) x y (M) ( 


hence, 


SPL 


N i 

1 m=0 


< r ) “ /(M) (r) 


{OM 

Ti^{ sp L \{f( R ^ s )} 


(27) 


iVl 
m =0 


(r) - / (re* 


v , 2tt 

x Y(k,£) I 


M 


By Proposition |T] we obtain the inequality: 


/ re*M s -5PL 


{/(P m e^)} 


iVi 

m=0 


(r) 


< C 


d 4 / (re^) 


dr 4 


x /i 4 . 


-^OO 


Hence, from (|27j) we obtain the inequality 

f ■ 

771—0 


SPL 


<7 C 

“ M 


< —C 
“ M 


= 2^C 


{/^ (p m )} 

d 4 / (re*) 


( r ) - /(M) ( r ) 


dr 4 

d 4 / (re lip ) 


M 


_ x/ > 4 E 

Loo(Dr) s=l 

M 


v l 27r 


M 


dr 4 
d 4 / (re**) 


x h 4 


^oo (Dll') 


7r 


dr 4 


x /i 4 . 


^oo (Dr) 


Finally, we obtain the final result: 

Tspline / r \ _ rpoly / r \ 

1 (N,M,K,N 1 ) ' (N,oo,K) > 


< 


I< a k N 


1 


fc=0f=lj=l 
K a k N 


2EEE^?-(M) x X {/(M) ^ m )} TO=0 (v^CM)) /(M) (v^h(M)) 

d 4 / (re^) 


<EEEE 

fc=0£=lj=l 


^h(M) x y ,(fc,t) 


x 2y/nC 


dr 4 


x h 4 . 


^oo (Dr) 
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By applying inequality (1221) we finish the proof. 


Remark 4 From the proof above we see that there is an alternative way to understand 
our scheme of obtaining the Hybrid Polyharmonic Cubature formula m we approximate 
the values of the function 


by means of the spline values SPL 


f ( \Ah(M) e 


{/ {jd-mP 1 M | (-\A?,(M)) > an d we use this ap¬ 
proximation to find • We see that numerically, the original scheme is 

more space-saving since we do not need to keep the approximations to the set of values 
{/ } - they are N X (2 K - 1) X M. 


In our experiments below we have chosen for simplicity 


R_ 

Ni' 


Rj = —j for j = 1, 2, Ni, 


(28) 


and we see that 


(/) - TPm, K) (/) < x ^ 


rP°ly 


H 4 


fir x 


d 4 f (re**) 


dr 4 


Loo {Dr) 


As proved in the following corollary, the Hybrid Polyharmonic Cubature formula 
is approximately exact for the same subspace of polynomials of the type r 2s r k Y/ k g\ (tp) 
which was considered in 


Corollary 5 For functions F = r 2s+k Y\^ (tp ), with 0 < s < 2 N — 1, k < M — 1 — K, and 
t = 1, and for Rj given by 128\) we have the estimate 


j-spline ( f) _ T (f) 


j^2s-\-k 

<cx 


where C is the constant provided by Proposition Q 
Remark 6 We see that for the weight w l 1 ) we have, m, 


ir (1) = - + cosp = ^—^(0,1) M + V^Y(i,i) (<f) 


hence, 


iv 


(i) 


r 1 r 1 2 

/ - rdr + / \Jirrdr = V2tt -\ — y/n ~ 3.4. 

Jo r J 0 2 


For the weight u/ 2 ) we obtain 


1 


(i) 2^/2 H) ( \ 4 

>4)( r ) = 7^ r ^(2M)W = -^4^T 


r for k > 1. 
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Hence 




2\/2 


2 dr + ^ 


k =l 



4 1 

4 k 2 — 1 


" 2 dr 


2y/2 


4 y> 1 < 2\/2 

3-^ 4k 2 — 1 ~ 30F ^ 

/c=l 




1 , 1 

—^- dx H— 

4x 2 — 1 3 


2\/2 4 

3-v/tF + 30F 


x 0.28 


0.742 


Corollary 0 shows that for R = 1, iVi ~ 10, and C ~ ^ (by Theorem 4.4.8 in 
[15]). for the weigths w^\ uJ 2 \ the formula I^'mk ivy (/) nearly exact with precision 
approximately given by 


C x 


R 2s+k „ „ 5 1 

-i— tw <-3.4 

Nf 11 11 ~ 384 10000 


4.4271 x 10 -6 . 


Hence, we may consider that the coefficient of efficiency of the cubature formula K (/) 
is approximately given by 

2AT x 2 (M - 1 - J\) 

~~ 3JVi x M ' 

We see that if we choose N\ ~ A r , and the Fourier approximation parameter K is given, 
then in order to make E k, 1 we have to choose 


M*tA(K + \) 

However, unlike the references [65], [66], we have to take care of the error bounds which 
depend on the derivatives of the functions / and the asymptotic properties of the weight 
functions hence we have to make even a more careful choice of the parameters M 

and K. 

3 Experimental results for the weight function (x, y ) 

As in [43]. we test the Hybrid Polyharmonic Cubature formula for integrals of the type 

/»27T pR 

I w (/) = / / / (r cos tp, r sin tp) ■ (r cos tp, r sin tp) ■ rdrdtp 

Jo Jo 

first for the weight function 

. 1 + x 1 \J2tt . . ,— , . 

' ( x , y) = 7 = = - + cos tp =- 1 ( 0 , 1 ) (tp) + V7rE (ljl) (tp). 

yr 2 + y r r 

We consider for our experiments three test functions: 
fo (x, y) = l+x 4 + y 3 

3 7 

fi(x,y) = 1 + ^ + T 2 j _„2 = 1 + r 2 cos 3 p -j- r 5 sin 7 tp 

i/r + y z x +y 

fi (x, y) = cos (ax + by) for a = 10 and b = 20. 

The first test function /o is a polynomial of degree 4, the second fi is not smooth at 0 
and the last one /2 is of oscillatory type. 
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3.1 Results for the Hybrid Polyharmonic Cubature Formula 

As in [43] . the number M of points on the circles is choosen to be equal to 9,25,63,83 
and the number N of concentric circles is chosen to be equal to 10,15, 25, 35, 50. The 
Gauss-Jacobi quadrature rules are with N\ = N points. We have the exact values given 
by 

h (/ow (1) ) = ^ti- « 6.754424205218060 
I] = —71- « 6.872 233 929727 67 

h (/ 2 w (1) ) « 0.301310 995335 215 


Below is the table with the values of K N] ) (/o), and the error. We see that the 

error is bigger compared with the similar table in [43] . 


N/M 

9 

25 

63 

83 

10 

6,754363639426710 

6,754363639426710 

6,754363639426710 

6,754363639426710 

15 6,754415757033810 

6,754415757033810 

6,754415757033810 

6,754415757033800 

25 

6,754423468244570 

6,754423468244570 

6,754423468244570 

6,754423468244570 

35 6,75424054924440 

6,754424054924440 

6,754424054924430 

6,754424054924440 

50 

6,754424177151970 

6,754424177151970 

6,754424177151970 

6,754424177151970 

Error 





N/M 

9 

25 

63 

83 

10 

-0,000060565791338 

-0,000060565791338 

-0,000060565791338 

-0,000060565791338 

15 

-0,000008448184238 

-0,000008448184238 

-0,000008448184238 

-0,000008448184248 

25 

-0,000000736973478 

-0,000000736973478 

-0,000000736973478 

-0,000000736973478 

35 

-0,000000150293608 

-0,000000150293608 

-0,000000150293618 

-0,000000150293608 

50 

-0,000000028066078 

-0,000000028066078 

-0,000000028066078 

-0,000000028066078 


Hybrid Polyharmonic cubature for fow w = (l + x 4 + y 3 ) (- + 


cos if 


43 

True value is —tt « 6.754424205218060 
20 


and the error. 


Below is the table with the values of k ) (■^ 1 ) 

is bigger compared with the similar table in [43] : 


and the error. Again, the error 


11 





















N/M 

9 

25 

63 

83 

10 

6,87224296287783 

6,87224296287783 

6,87224296287783 

6,87224296287783 

15 

6,87223588060173 

6,87223588060173 

6,87223588060173 

6,87223588060173 

25 

6,87223420205342 

6,87223420205342 

6,87223420205342 

6,87223420205342 

35 

6,87223400297000 

6,87223400297000 

6,87223400297000 

6,87223400297000 

50 

6,87223394775545 

6,87223394775545 

6,87223394775545 

6,87223394775545 

Error 





N/M 

9 

25 

63 

83 

10 

0,00000903315016 

0,00000903315016 

0,00000903315016 

0,00000903315016 

15 

0,00000195087406 

0,00000195087406 

0,00000195087406 

0,00000195087406 

25 

0,00000027232575 

0,00000027232575 

0,00000027232575 

0,00000027232575 

35 

0,00000007324233 

0,00000007324233 

0,00000007324233 

0,00000007324233 

50 

0,00000001802778 

,00000001802778 

0,00000001802778 

0,00000001802778 


Hybrid Polyharmonic cubature for f \'«/ 1 - 1 = (l + r 2 cos 3 p + r 5 sin' — f cos p^ 


35 

True value is — it ~ 6. 872 233 929 727 67 
16 


There is no mistake in the table above since results do not change with Ml This is 
due to the fact that 


■ 7 7 21 1 . 7 , 35 . 

sm P = -T77 Sin bp - — sin 3 p - — sin 7 p+ — sin p 

64 64 64 64 

which implies that the function g = /i?i/b has no terms g^/) ( r ) with k > 4; hence, the 
interpolation cubic splines of the polynomials g^k,e) (r) f coincide with them. 

For the oscillatory test function we obtain 


N\M 

9 

25 

63 

83 

10 

0,19210087475239 

0,58134400368821 

0,56846433865624 

0,56846433865624 

15 

-0,00842490123728 

0,38518390970894 

0,37239275649383 

0,37239275649383 

25 

-0,08069670830424 

0,31431334300881 

0,30152604401835 

0,30152604401835 

35 

-0,08150067781664 

0,31360790761542 

0,30081999553130 

0,30081999553130 

50 

-0,08116599113863 

0,31395572503057 

0,30116759220177 

0,30116759220177 


Error 


N\M 

9 

25 

63 

83 

10 

-0,10921012058282 

0,28003300835300 

0,26715334332102 

0,26715334332102 

15 

-0,30973589657249 

0,08387291437372 

0,07108176115862 

0,07108176115862 

25 

-0,38200770363946 

0,01300234767360 

0,00021504868313 

0,00021504868313 

35 

-0,38281167315186 

0,01229691228021 

-0,00049099980392 

-0,00049099980392 

50 

-0,38247698647384 

0,01264472969535 

-0,00014340313345 

-0,00014340313345 


Hybrid Polyharmonic cubature for f2W^ = (cos (lOx + 20y)) ( —I- cos p 


True value is h (/2W (1) ) ~ 0.301310 995 335 215 
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4 Experimental results for the weight function w ^ (x, y ) 

As in |43], we consider the second weight function: 

u/ 2 ) ( re lv ) := |y| = |rsiny>| 

We have its orthonormalized Fourier coefficients: 


2^/2 ^ ^ 

^(0,1) (0 = and ™(2M) (0 = ~^ 4 k 2 _ I r for ^ L 

5 Results for the Hybrid Discrete Polyharmonic Cubature 

As in |43j, we consider the test functions 

h {x,y) = 30x 12 , h(x,y) = \y \, 


(29) 


for which 


4(2) (/) = Ii (30x 12 |y|) = « 0.615 384615 384 616 

-L O 

n 27T ^ 

|r sin 4 2 diprdr = -tt « 0.785 398 163 397 448 

The table contains the results for the Hybrid Discrete Polyharmonic Cubature 
for the test function 4 (x,y) (fMl) . where Ji = 22: 


N\M 

9 

25 

63 

83 

10 0,565617343585166 

0,620572422003199 

0,620572422003199 

0,620572422003199 

15 

0,561709413462587 

0,616243839415133 

0,616243839415133 

0,616243839415132 

25 0,561006480042405 

0,615463257008362 

0,615463257008361 

0,615463257008362 

35 

0,560949764331337 

0,615400379071044 

0,615400379071044 

0,615400379071044 

50 

0,560937835197964 

0,615387283068315 

0,615387283068315 

0,615387283068316 


Error 


N\M 

9 

25 

63 

83 

10 

-0,0497672717994500 

0,0051878066185830 

0,0051878066185830 

0,0051878066185830 

15 

-0,0536752019220290 

0,0008592240305171 

0,0008592240305171 

0,0008592240305161 

25 

-0,0543781353422109 

0,0000786416237460 

0,0000786416237450 

0,0000786416237460 

35 

-0,0544348510532789 

0,0000157636864280 

0,0000157636864280 

0,0000157636864280 

50 

-0,0544467801866519 

0,0000026676836991 

0,0000026676836991 

0,0000026676837001 


Hybrid Discrete Polyharmonic cubature for 30x 12 \y\ 


True value is m 0.615 384 615 384616 

For the test function 4 ( x,y ) = \y \, with IF = 22 in the formula (1241) . we obtain the 
table 
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N\M 

9 

25 

63 

83 

10 

0,785206660 

0,785352337 

0,785367124 

0,785369362 

15 

0,785208297 

0,78535897 

0,785373081 

0,785375274 

25 

0,785208235 

0,785361119 

0,785374994 

0,785377171 

35 

0,785208149 

0,78536144 

0,785375276 

0,785377452 

50 

0,785208109 

0,785361541 

0,785375364 

0,785377539 


Error 


N\M 

9 

25 

63 

83 

10 

-0,0001915037 

-0,0000458267 

-0,0000310393 

-0,0000288009 

15 

-0,0001898666 

-0,0000391939 

-0,0000250824 

-0,0000228890 

25 

-0,0001899281 

-0,0000370443 

-0,0000231697 

-0,0000209919 

35 

-0,0001900142 

-0,0000367231 

-0,0000228870 

-0,0000207116 

50 

-0,0001900544 

-0,0000366227 

-0,0000227993 

-0,0000206248 


Hybrid Polyharmonic cubature for = |y| 2 


True value is as 0.785 398 163 397448 


6 Conclusions 

The comparison with the results in [35] shows that the Hybrid Polyharmonic Cubature For¬ 
mula provides a satisfactory tradeoff between precision and number of evaluation points. 
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